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Disconnect ivity graphs are used to characterize the potential energy surfaces of Lennard-Jones 
clusters containing 13, 19, 31, 38, 55 and 75 atoms. This set includes members which exhibit either 
one or two 'funnels' whose low-energy regions may be dominated by a single deep minimum or contain 
a number of competing structures. The graphs evolve in size due to these specific size effects and 
an exponential increase in the number of local minima with the number of atoms. To combat the 
vast number of minima we investigate the use of monotonic sequence basins as the fundamental 
topographical unit. Finally, we examine disconnectivity graphs for a transformed energy landscape 
to explain why the transformation provides a useful approach to the global optimization problem. 



I. INTRODUCTION 

The relationship between dynamics and the potential 
energy surface (PES), or potential energy 'landscape', 
holds the key to understanding a wide range of phenom- 
ena, including protein folding, global optimization and 
the properties of glasses. For example, what are the to- 
pographical features of the PES that enable a protein to 
fold to its native state? What sort of PES renders global 
optimization tractable, and how can it be transformed!!! 
to make global optimization easier? Which features of 
the PES produce good glass formers, and what distin- 
guishes 'strong' from 'fragile' liquids^? 

Stillinger and Webem first showed the utility of parti- 
tioning configuration space into basins of attraction sur- 
rounding each local minimum. The thermodynamics of 
the system can then be forpjuilated in terms of the proper- 
ties of these local minima.QQ However, the minima alone 
contain no dynamic information. Therefore, to obtain in- 
sight into the dynamics it is necessary to locate the tran- 
sition states which connect the local minima. This infor- 
mation can then be incorporated into a master equation 
that describes the flow of probability between minima 
as the system evolves towards equilibrium.ETEl As well 
as providing a picture of the relationship between the 
dynamics and the PES, this approach can probe time 
scales muchJonger than those accessible in conventional 
simulations .□ 

We also wish to relate the thermodynamics and dy- 
namics to general features of the PES in a qualitative 
manner. Disconnectivity graphs are a recently-developed 
tool that provide a helpful visualization of the PES. 
These graphs..haiie now been calculatecLfop-a number of 
polypeptidesolld'EJ and various clusters .E3~EJ They show 
which of the minima in a sample arc connected by path- 
ways lying below a series of energy thresholds, and so 
provide a picture of the hierarchy of energy barriers in a 
system. 

In this paper we present disconnectivity graphs for a 
series of Lennard-Jones (LJ) clusters. Our aims are to 
gain more insight into the potential energy landscapes 
of these clusters and to explore systematically how the 
graphs evolve with size. LJ clusters provide an ideal test 



case because their dynamics and thermodynamics have 
been thoroughly studied and they have been usedria the 
development of global optimization algorithms .E-TLJ 

As the size of the cluster increases the number of min- 
ima and transition states is expected to grow at least 
cxponentialhz- even when permutational isomers are not 
included£3~E3 Furthermore, cluster properties do not 
vary smoothly with the number of atoms in this size 
regime, and we have chosen to examine sizes which should 
illustrate particularly interesting features. 

Due to the exponential increase in the number of sta- 
tionary points with size it is difficult for a disconnec- 
tivity graph to retain a global picture of the PES. We 
have therefore explored the use of monotonic sequence 
basinscil as the funda menta l topographical unit r ather 
than minima (Section III B ) . Finally, in Section III C 



we use disconnectivity graphs to probe the transformed 
energy landscape that is searched by the Monte Carlo 
minimizationE3 or basin-hoppingtH algorithm. 



II. METHODS 



The atoms interact via a Lennard-Jones potential 
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where e is the pair well depth and 2 1 / 6 ct is the equilibrium 
pair separation. 

To examine the topography of the PES, we need to 
locate its minima and the network of transition states 
and pathways that connect them. Transition sts 
be found efficiently using eigenvector-following! 
which the energy is maximized along one direction and 
simultaneously minimized in all the others. The min- 
ima connected to a transition state are defined by the 
two steepest-descent paths which begin parallel and an- 
tiparallel to the unique Hessian eigenvector whose cor- 
responding eigenvalue is negative. The steepest-descent 
paths were calculated usisg a method which employs an- 
alytic second derivatives. E3 
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The samples of minima and transition states were gen- 
erated pHsie,s-an approach which has also been described 
before Jl3 ! tiHla The elementary step in this process is to 
perform a transition state search after stepping away 
from a minimum along one of the Hessian eigenvectors 
and then, if this search finds a transition state, to gen- 
erate the corresponding pathway. If the transition state 
is connected to the minimum from which the search was 
started, we add it to our database along with the con- 
nected minimum. Occasionally we find that the transi- 
tion state is not connected to the starting minimum. In 
this case we only add it to the database if it is connected 
to a minimum that is already present. At each step we 
start from the lowest-energy minimum from which fewer 
than a specified number, 2n ov , of transition state searches 
have been performed. Searches are carried out in positive 
and negative directions along each eigenvector in order 
of increasing eigenvalue. The value chosen for n cv deter- 
mines how thoroughly the PES is searched in the vicinity 
of a given minimum. 

By repeating this process until no new minima are 
found we can obtain a nearly exhaustive catalogue of 
the minima. This approach was used for the 13-atom 
cluster, LJ13, but for the larger clusters it is impossible 
to find all the minima. We therefore stopped searching 
once we were confident that we had obtained an accu- 
rate representation of the low-energy regions of the PES. 
For LJ75 an alternative strategy was required due to the 
double-funnel character of the landscape; this approach 
is described in the following section. 




FIG. 1. A selection of low-energy minima for LJ13, LJ19, 
LJ31, LJ38, LJ55 and LJ75. The first number denotes the 
number of atoms, and the second denotes the energetic rank 
of the minimum (the global minimum has rank 1, etc.). 



III. RESULTS 

We now characterize the PES's of LJ 13 , LJ 19 , LJ 31 , 
LJ38, LJ55 and LJ75. To explain this selection of clus- 
ter sizes, and as a background to the interpretation of 



the disconnectivity graphs, we briefly review some of the 
structural properties of small LJ clusters. Comparisons 
of particularly stable sequences of clusters indicate that 
structures based on Mackay icosahedraES (e.g.riS.l and 
55.1 in Figure are dominant up to N ~ 1600.E3 There- 
fore, as can be seen from Figure |2j, most of the clusters in 
the size range that we are considering have a global min- 
imum based on icosahedral packing. The sizes at which 
complete Mackay icosahedra are geometrically possible 
(N=13, 55,. . .) are particularly stable compared to other 
sizes, leading to particularly low-energy global minima 
(Figure ||). For these clusters, there is a large energy gap 
between the two lowest-energy minima (2.85e and 2.64e 
for LJ13 and LJ 55 , respectively) and we expect the PES 
to have a single deep funnel. 

The term 'funnel' was introduced in the protein fold- 
ing literature to describe the situation where a collection 
of downhill pathways converge on the native structure 
of the protein .EH As a result of this PES topography 
the protein is 'funnelled' towards the native structure on 
relaxation. We use the term here in a similar manner 
except that, instead of converging on the native struc- 
ture of the protein, the funnel should converge on a sin- 
gle low-energy minimum, or possibly to a collection of 
structurally-related low-energy minima. 




FIG. 2. Comparison of the energies of icosahedral (I), 
decahedral (D) and close-packed (C) LJjv clusters. The en- 
ergy zero is -Emi, a function fitted to the energies of the 
first four Mackay icosahedra at JV=13, 55, 147 and 309: 
E M i = -2.3476 - 5.4633iV 1/3 + 14.8814iV 2/3 -8.5699JV. The 
sizes that we consider in detail in this paper have been la- 
belled. 

Between N—13 and 55 the energy first increases rel- 
ative to an energy function interpolated between the 
Mackay icosahedra as an overlayer is added to the 13- 
atom icosahedron, and then decreases as the overlayer 
approaches completion at iV=55 (Figure ||). Of course, 
there are smaller features superimposed on this broad 
maximum. At certain sizes particularly stable overlayers 
with no low-coordinate atoms can be formed, giving rise 
to subsidiary minima in Figure^. For example, at N=19 
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a six-atom cap on the 13-atom icosahedron leads to the 
relatively stable double icosahedron shown in Figure [l]. 
We would expect LJig to exhibit a single-funnel PES but 
to relax less efficiently to the global minimum than LJ13 
or LJ55. 

We also chose to study a size, N=31, from near the 
top of the broad maximum in Figure |^, where there are a 
large number of low-energy minima associated with the 
various ways that the atoms of the overlayer on the 13- 
atom icosahedron can be arranged. Moreover, for this 
size the lowest-energy decahedral cluster is only a little 
higher in energy than the icosahedral structures. 

For most sizes the global minimum has icosahedral 
character, but there are a few clusters for which this 
is not the case. For N=38 the global minirniun_is the 
face-centred-cubic (fee) truncated octahedrontSE3 (38.1 
in Figure [l]) and for 75 < N < 77 .the global minima 
are based on the Marks decahedronEHl (75.1 in Figure 
|l|). We chose to study iV=38 and 75 where the global 
minimum is structurally very different from the second 
lowest-energy minimum, which in each case is icosahe- 
dral. The two lowest-energy minima are therefore rather 
distant in configuration space, and they are_separated 
by large potential energyca and free energyo barriers. 
Therefore, we expect these surfaces to have two funnels: 
one which leads to the low-energy icosahedral minima 
and one which leads to the global minimum. 

Details of the samples of minima and the transition 
states that we found for each of the clusters are given in 
Table Q. For LJ13, we were able to gauge how close the 
samples were to convergence by following the number of 
minima and transition states as n cv was increased. The 
number of minima seemed close to reaching an asymp- 
tote, but the number of transition states was still increas- 
ing at n cv =15.Q 



A. Disconnectivity graphs 

The conceptual basis for the disconnectivity graph ap- 
proach is as follows.EEl At a given total energy, E, min- 
ima can be grouped into disjoint sets, called superbasins, 
whose members are mutually accessible at that energy. 
In other words, each pair of minima in a superbasin are 
connected directly or through other minima by a path 
whose energy never exceeds E 7 but would require more 
energy to reach a minimum in another superbasin. At low 
energy there is just one superbasin — that containing the 
global minimum. At successively higher energies, more 
superbasins come into play as new minima are reached. 
At still higher energies, the superbasins coalesce as higher 
barriers are overcome, until finally there is just one con- 
taining all the minima (provided there are no infinite 
barriers). 



TABLE I. The number of minima n m i n , transition states 
n ts , and monotonic sequence basins iimsb in our samples. 
For each sample, transition state searches were performed 
from the lowest n soarc h minima, and for these minima tran- 
sition state searches were performed along the lowest n ev 
eigenvalues. Values of timsb and n SGarc h are not included for 
LJ75 because a different method was used to search the PES. 
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To construct a disconnectivity graph, the superbasin 
analysis is performed at a series of energies. At each 
energy, a superbasin is represented by a node, and lines 
join nodes in one level to their daughter nodes in the level 
below. The horizontal position of a node has no signifi- 
cance, and is chosen for clarity. Every line terminates at 
a local minimum. 

The disconnectivity graphs produced using this proce- 
dure are shown in Figure ||. For A=13 it is possible to 
include all the minima that we have found in the discon- 
nectivity graph, which therefore provides a practically 
complete global picture of the PES. The graph has the 
form expected for an ideal funnel: there is a single stem, 
representing the superbasin of the global minimum, with 
branches sprouting directly from it at each level, indicat- 
ing the progressive exclusion of minima as the energy is 
decreased. 

The form of this graph implies that all the minima 
are directly connected to the superbasin associated with 
the global minimum. In fact, 99% of the minima are 
within two rearrangements of the global minimum and 
none are further than three steps away. Furthermore, 
there are 911 structurally distinct transition states con- 
necting 535 minima directly to the global minimum. If 
the reaction path degeneracy is taken into account, there 
are 108 967 transition states (some of which are permuta- 
tional isomers) connected to each permutational isomer 
of the global minimum.L3 

If we examine the distance of a minimum from the 
global minimum as a function of its potential energy, for 
LJ13 the resulting plot (Figure ^a) has the form we would 
expect of a funnel: the distance from the global mini- 
mum increases as the potential energy of the minimum 
increases. It is also worth noting that the slope of the fun- 
nel is steeper than for any of the other clusters. Hence it 
is no surprise that relaxation down the PES to the global 
minimum is relatively easy for this cluster. However, re- 
laxation is hindered somewhat by the fairly large barriers 
(l-2e) that exist for escaping from some of the minima 
into the superbasin of the global minimum (Figure ^a). 
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FIG. 3. Disconnectivity graphs for (a) LJ13, (b) LJ19, (c) LJ31, (d) LJ38, (e) LJ55 and (f) LJ75. In (a) all the minima 
are represented. In the rest only the branches leading to the (b) 250, (c) 200, (d) 150, (e) 900 and (f) 250 lowest-energy 
minima are shown. The numbers adjacent to the nodes indicate the number of minima the nodes represent. The branches 
associated with the minima depicted in Figure ^ are labelled by their energetic rank. In (e) an enlarged view of the branch 
marked i is shown in the inset. The energy scale is in units of e. 
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FIG. 3. cont 



For all the other clusters it is not computationally fea- 
sible to obtain a nearly complete set of minima. More- 
over, if we attempt to represent all the minima of our 
samples on a disconnectivity graph, the density of lines 
simply becomes too great. Instead we only represent 
those branches that lead to a specified number of the 
lowest-energy minima, both for clarity and because our 
samples of minima and transition states are likely to be 
most complete for the low-energy regions of the PES. We 
should note that the minima that are not represented 
can still contribute to the appearance of the graph if 
they mediate low barrier paths between minima that are 
included. However, this approach does have the conse- 
quence that, as the size increases, the graphs increasingly 
focus on a smaller and smaller proportion of the surface. 

We can see this effect if we examine the graph of L J55 
for which we also expect the PES to exhibit a single fun- 
nel. Although there is more fine structure, the form is 
similar to that of LJ13. However, all the minima rep- 
resented in the disconnectivity graph have relatively or- 
dered structures and so, unlike the graph for LJ13, this 
representation does not tell us whether there is a fun- 
nel leading from the disordered liquid-like minima to 
the global minimum Mackay icosahedron. Instead, it 
only shows that the low-energy region of the PES asso- 
ciated with structures based on the Mackay icosahedron 
is funnel-like. The graph probably only represents the 



bottom of a larger funnel leading down from the huge 
number of disordered minima. The plot of the distance 
from the global minimum as a function of potential en- 
ergy shows a glimpse of this larger funnel (Figure |^e) . 

It is easy to understand the differences between the 
graphs for LJ13 and LJ55. The energy gap between the 
global minimum and the disordered liquid-like minima 
is an extensive quantity, whereas the energy gap be- 
tween two ordered structures, where a few atoms have 
changed position, is related to the change in the num- 
ber of nearest-neighbour contacts and is independent of 
size. Furthermore, the number of possible ways of ar- 
ranging surface atoms and vacancies in an ordered struc- 
ture becomes increasingly large. Therefore, the number 
of ordered minima which lie below the liquid-like band of 
minima increases rapidly with the number of atoms. 

The fine structure of the LJ55 disconnectivity graph 
reveals some interesting features. The minima separate 
into bands related to th|&,number of defects present in 
the Mackay icosahedron£3 For example, all the minima 
in the first band above the global minimum are Mackay 
icosahedra with a missing vertex and an atom on the 
surface. The eleven minima in this band correspond to 
the eleven possible sites for this atom that are unrelated 
by symmetry Of these minima the four lowest in en- 
ergy have the atom located over the centre of one of the 
faces of the Mackay icosahedron, and the other seven 
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have the atom off-centre. In the disconnectivity graph 
these minima split into four groups corresponding to the 
four symmetry-unrelated faces on which the adatom can 
be located. The splitting occurs because the barriers for 
rearrangements in which the adatom passes between faces 
are larger than the barriers for changing the position of 
the adatom on a face. Therefore, in the disconnectivity 
graph the minima first become connected to the other 
(one or two) minima with the adatom on the same face 
before they are connected to the stem associated with 
the global minimum. 

The second band of minima consists of Mackay icosahe- 
dra with two missing vertices and two surface atoms. The 
lower-lying minima in this band have the two adatoms 
in contact, either on the same face or bridging an edge. 
The minima with the two adatoms unpaired give rise to 
a repeated motif, an example of which is illustrated in 
the inset. This feature is repeated over ten times on the 
left-hand side of the graph with the top node always at 
— 273. Oe. The minima split into these sets because of the 
lower barriers for an adatom moving between sites on the 
same face. If the two faces with adatoms are unrelated 
by symmetry, there are sixteen distinct ways of arranging 
the atoms on the two faces. The lowest-energy minimum 
of these sixteen has the two adatoms in the central sites. 
Slightly higher in energy are the six minima with one 
adatom central and one off-centre. Higher still are the 
nine minima with both adatoms off-centre. 

The features noted above are reflected in the dynamics 
of L J 55 . Just below the melting point the cluster passes 
between the Mackay icosahedron and states with one and 
two defects, residing in each state for periods of the order 
of nanoseconds (using potential parameters appropriate 
to argon) While in one of these states, the surface 
atoms and vacancies diffuse across the surface. This time 
scale separation between diffusive motion and the forma- 
tion and annihilation of defects is a consequence of the 
higher barriers for the latter processes and the distances 
over which a surface atom and vacancy must diffuse to 
recombine. The latter is reflected in the wide range of 
distances that minima in this second band are away from 
the global minima. 

The two minima (55.13 and 55.149) in Figure ||e that 
must overcome the highest barriers to reach the super- 
basin containing the global minimum do not have icosa- 
hedral structures. Both have D$h point group symmetry 
and are depicted in Figure U The lower -energy of the two 
is actually the thirteenth lowest-energy minimum and 
was recently discovered by Wolf and Landman. c3 The. 
other minimum has previously been noted by Wales SB 
Both can be converted into a Mackay icosahedron by a 
single cooperative rearrangement in which parts of the 
structure twist around the fivefold axisO As a signifi- 
cant number of nearest-neighbour contacts are broken in 
these rearrangements, the barriers are higher than those 
for the localized rearrangements by which the defective 
Mackay icosahedra interconvert. 

The disconnectivity graph for LJ19 shows that the 



PES is again funnel- like (Figure |^b), as expected. The 
branches for most of the minima connect directly to the 
superbasin containing the global minimum. Although 
our graph only shows branches leading to the lowest 250 
minima, Figure ^b reveals that the funnel continues up 
to higher energies. 

The PES jOf_LJj g has been analysed previously in sev- 
eral studies J3 : E1l3 The profiles of the downhill pathways 
to the global minimum were described as 'sawtooth-like' 
rather than 'staircase-like' because the barrier heights are 
relatively large compared to the energy differences be- 
tween the minima. On this basis Ball et al. concluded 
that LJ19 has topographical features typical of a glass- 
former. The disconnectivity graph also shows that some 
of the downhill barriers are quite large. However, as the 
global minimum is at the bottom of a funnel the barriers 
only slow down the rate of relaxation towards the dou- 
ble icosahedron, rather than preventing it. Lowering the 
energy does not take the system away from the global 
minimum, but rather towards it (Figure ^b). The ap- 
pearance of N = 19-as a magic number in mass spectra 
of rare gas clusterscil confirms that the global minimum 
is kinetically accessible. 

The disconnectivity graph for LJ31 (Figure |3|c) is fun- 
damentally different from those we have considered so 
far. The energetic bias towards the global minimum is 
smaller (Figure^,). In fact there are a number of min- 
ima with energies close to that of the global minimum 
which are separated from it by fairly large energy bar- 
riers. This situation is partly the result of competition 
between two distinct types of overlayer. In the first type, 
the anti-Mackay overlayer, atoms are added to the faces 
and vertices of the underlying 13-atom icosahedron (giv- 
ing rise, for example, to the double icosahedron, 19.1). In 
the second type, the Mackay overlayer atoms are added 
to the edges and vertices. The completion of the Mackay 
overlayer leads to the next Mackay icosahedron. LJ31 is 
the first size for which a cluster with the Mackay over- 
layer is the global minimum. c3 It can be seen from Figure 
gthat minimum 31.1 is a fragment of the 55-atom Mackay 
icosahedron. The second lowest-energy minimum, 31.2, 
has an anti-Mackay overlayer. 

There are also some low-energy decahedral minima for 
LJ31. For example, minimum 31.9 (Figure [l]) is separated 
by a large energy barrier from the global minimum, be- 
cause not only must there be a change in morphology 
from decahedral to icosahedral, but the cluster must also 
change shape. Some of the decahedra with more spher- 
ical shapes are connected to the superbasin associated 
with the global minimum by smaller barriers. 

We can deduce something of the relaxation behaviour 
of LJ31 from its disconnectivity graph. Once the clus- 
ter has reached a low-energy configuration, presumably 
by rapid descent of a funnel from the liquid, subsequent 
relaxation towards the global minimum may be consid- 
erably slower. There is little energetic bias at the bot- 
tom of the PES to guide the system towards the global 
minimum and the barriers for intcrconversion of the low- 
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energy minima can be relatively large. Therefore, it is 
not surprising that the time required to find the global 



minimum using the basin-hopping global optimization al- 
gorithm shows a maximum at iV=31.E3 




The effects of competing structures that we noted 
for LJai appear in a more extreme form for LJ38 and 
L J75 .M For both these clusters the disconnectivity graph 
clearly separates the low-energy minima into two main 
groups, namely those associated with the global mini- 



mum and those with icosahedral structure. These two 
groups of minima are separated by a large energy bar- 
rier, so the graph splits into two stems at high energy 
which lead down to two structurally distinct sets of low- 
energy minima. This splitting is characteristic of a mul- 
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tiple funnel PES. The separation is particularly dramatic 
for LJ75, where the decahedral to icosahedral barrier is 
over 3e larger than any of the other barriers between the 
250 lowest-energy minima. Although the disconnectivity 
graphs clearly show only the bottom of the PES, these 
two groups of minima can be associated with separate 
funnels, and give rise to distinct thermodynamic states.E-3 
The double funnel structure is also apparent from the 
plots in Figure [l|d and f. Passing from the global mini- 
mum to the icosahedral funnel, the energy first increases 
as the primary funnel is ascended and then decreases dur- 
ing the descent into the second funnel. 

From the relative numbers of minima associated with 
each funnel it is clear that the icosahedral funnel is much 
wider. This is one of the reasons why relaxation is much 
more likely to lead to the icosahedral minima. There are 
also thermodynamic effects: the region containing the 
icosahedral minima becomes lowestin-iree energy at low 
temperature (T^O.Uek- 1 for LJasBEii and T^Qmek- 1 
for LJ75EZI) because of the entropy that arises from the 
large number of low-energy icosahedral minima. Between 
this transition temperature and the melting point there 
is therefore a thermodynamic driving force towards the 
icosahedral structures. Moreover, for LJ38 we have shown 
that the free energy barrier for entering the icosahedral 
region of configuration space is lower than for entering 
the fee funnel.li-3 

Once the cluster enters the icosahedral funnel it is 
likely to be trapped because of the large energy and free 
energy barriers to passing between the two funnels. The 
energy barriers are 4.22e and 3.54e for LJ 38 and 8.69e and 
7.48e for LJ75. At higher temperatures the entropy of the 
intermediate states reduces theJree energy barriers from 
these zero temperature limits.Efl For the above reasons 
LJ38 and LJ75 provide relatively difficult test cases for 
any unbiased global optimization algorithm. It is only 
relatively recently that algorit hms, have begun to find 
the truncated pctahedron,Ej1i3't§E3~E3 and only one utu. 
biased methodO and one method that involves seedingEa 
have reported finding the Marks decahedron. 

The difficulty in finding the global minimum for these 
clusters is illustrated by statistics for the basin- hopping 
algorithmO The time required to find the global min- 
imum for LJ 38 is a maximum with respect to neigh- 
bouring sizes and for LJ7 5 _77 the time is so long that 
good statistics for the first passage time have not yet 
been obtained. Global optimization is an order of mag- 
nitude more difficult for LJ75 than LJ38 for a combina- 
tion of reasons. First, LJ75 has a much larger search 
space and a much greater number of minima than LJ38. 
Second, the temperature at which the global free energy 
minimum becomes associated with the global potential 
energy minimum lies further below the melting point 
(T m «0. 17 efc- 1 for LJ38H and I^wO^efe- 1 for LJ 75 ). 
Third, the path between the global minimum and the 
lowest-energy icosahedral structure is also longer, more 
complicated and higher in energy than for LJ38. The 
transition involves not only a change in morphology, but 



also a change in shape — the Marks decahedron is oblate 
whereas the icosahedral structures are prolate. The path- 
way therefore involves both cooperative rearrangements, 
where the structure twists around a quasi- fivefold axis, 
and surface diffusion steps. The lowest-energy pathway 
that we found between the two lowest-energy minima is 
depicted in Figure |a. It passes through sixty-five tran- 
sition states. There are significantly shorter paths be- 
tween the same two minima (the shortest is 41. 5c), but 
these all involve higher barriers. There is also a signifi- 
cant difference in the character of the pathways between 
the minima at the bottom of the two funnels for LJ38 
and L J75 . All the structures along the L J75 pathway are 
ordered and solid-like, whereas at its highest points ;|| lhe 
LJ38 pathway passes through disordered structures IH 



For LJ75 the highest-energy minimum on the lowest- 
energy interfunnel path lies at position 10 909 in terms 
of an energy ranking of the minima in our sample, where 
1 is the global minimum. Many more low-lying minima 
could have been found if we had searched the icosahedral 
region of configuration space more intensively. There- 
fore, if we had simply performed consecutive transition 
state searches from minima in order of their energetic 
rank it is unlikely that a pathway between the two fun- 
nels could have been found — rather the search would have 
most likely been stuck in one of the funnels cataloguing 
the multitude of ordered structures with that morphol- 
ogy. To find a pathway connecting the two lowest-energy 
minima we therefore had to bias the search to probe inter- 
mediate structures. We first constructed a series of deca- 
hedral minima which were decreasingly prolate and then 
increasingly oblate. We then started searching the PES 
around these structures with the aim of finding pathways 
connecting them both to the Marks decahedron and to 
the low energy icosahedral structures. 



Once we had found a pathway between the two struc- 
tures, we only performed transition state searches from 
those minima that were connected to either of the two 
lowest-energy minima by a pathway that had no tran- 
sition states higher than the highest-energy transition 
state on the lowest-energy path between structures 75.1 
and 75.2. Moreover, this search concentrated on those 
minima in the set that had-aa intermediate value of the 
bond-order parameter, QgEaEj By this procedure, in- 
creasingly low-energy paths were found between the two 
funnels. There is, of course, no guarantee that we have 
found the lowest-energy pathway, but we doubt if a sig- 
nificantly lower one exists. 
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FIG. 5. Properties of the LJ75 PES. (a) The low- 
est-barrier path between the global minimum and the second 
lowest minimum. Each maximum corresponds to a transi- 
tion state, (b) Scatter plot of the bond order parameter, Qe, 
versus potential energy for all the minima in our sample, (c) 
Scatter plot of Qe for each minimum versus the energy of the 
highest transition state on the lowest-energy path to either 
structure 75.1 or 75.2, whichever path is the lower. The two 
+ symbols indicate the values of Q§ and potential energy for 
the two lowest-energy minima. 

From Figure |^b we can see how the bond-order param- 
eter can separate the minima into various groups. The 



group with Qg«0.02 correspond to icosahedral minima 
with an anti-Mackay overlayer (such as 75.2), whereas 
Q6=0.306 for the Marks decahedron. The decahedral 
minima based upon 75.1 have similar values, these be- 
ing generally lower for structures that are less oblate. 
There are also bands of minima at intermediate val- 
ues of Qq. For example, icosahedral structures with a 
Mackay overlayer have Qq~0.15. Some of the other min- 
ima with intermediate values have in part motifs sim- 
ilar to 55.13, and are connected to icosahedral or dec- 
ahedral minima by rearrangements involving concerted 
twists around quasi-fivefold axes. 

This distribution of Qq values provides us with another 
way of visualizing the double-funnel structure of the L J75 
PES. Figure ||c maps out the energy of the transition 
state at the top of the lowest-barrier pathway from each 
minimum to either 75.1 or 75.2, whichever is lower, as a 
function of Q e for the minimum. The plot shows that 
the barrier separating a decahedral minimum from the 
Marks decahedron increases as the value of Qq deviates 
further from the value for the global minimum. This 
trend continues until Q6~0.16. Below this value some 
of the minima have lower barriers for paths to minimum 
75.2, and the energies generally decrease as the value of 
Qe approaches that of 75.2 from above. The icosahedral 
minima with a Mackay overlayer stand out in the plot as 
they have both an intermediate value of Q e and low bar- 
riers connecting them to minimum 75.2. The transforma- 
tions between the two types of icosahedral overlayer are 
usually achieved by a concerted twisting of the overlayer 
around one of the vertices of the underlying icosahedron. 



B. Coarse-graining the PES 

As the size of the cluster increases our disconnectiv- 
ity graphs focus on an increasingly small proportion of 
the whole PES to avoid being swamped by the rapidly 
increasing number of minima. However, it would be de- 
sirable to retain a more global picture of the PES. To do 
so, the disconnectivity graphs would need to be based not 
on the barriers between minima, but between larger topo- 
graphical units. For example, Levy and Becker 'diluted' 
their sample of hexapeptide conformations by removing 
conformations that were energetically and structurally 
similar ILil However, this approach requires a meaningful 
measure of similarity, which is probably harder to devise 
for a cluster than a molecule with a bonded framework. 

Here, we explore the use of monotonic seq uen ces, a 
concept introduced by Berry and coworkers ,LfEilc3'Eil to 
produce a more coarse-grained picture of the PES. Mono- 
tonic sequences are series of connected minima where the 
potential energy decreases with every step. The collec- 
tion of sequences leading to a particular minimum de- 
fines a 'basin'. To avoid confusion with the various other 
'basins' that have been defined, we will always refer to 
such a set as a monotonic sequence basin (MSB). The 







MSB leading to the global minimum is termed 'primary', 
and is separated from neighbouring 'secondary' MSB's 
by transition states lying on a 'primary divide', and so 
on. It is important to realise that, above such a divide, 
it is possible for a minimum to belong to more than one 
MSB through different monotonic sequences. 

For the division of the PES provided by a monotonic 
sequence analysis to be useful in an analysis of the dy- 
namics, transitions between minima in an MSB must 
be more rapid than transitions between different MSB's. 
Kunz and BertivJpund evidence for this in a simplified 
model of LJigBEil However, there is no guarantee that 



this separation will always hold, because MSB's are de- 
fined by the connections between minima without refer- 
ence to the size of the barriers between them. 

Disconnectivity graphs that only include the minima 
at the bottom of each MSB can be produced by exclud- 
ing all the minima directly connected to a lower-energy 
minimum. In the resulting graphs, branches are joined 
by a node at the energy of the lowest transition state on 
the divide between the MSB's. For the clusters we con- 
sider here, the number of MSB's is small enough (Table 

for all of them to be represented on the graph. Two 
examples are shown in Figure ^. 



-167.5 , (a) 



-168.0 
-168.5 
-169.0 
-169.5 
-170.0 
-170.5 
-171.0 
-171.5 
-172.0 
-172.5 
-173.0 
-173.5 
-174.0 




280 



FIG. 6. Disconnectivity graphs for (a) LJ38 and (b) LJ55. Only the branches corresponding to monotonic sequence basin 
bottoms are shown. The branches associated with the minima depicted in Figure |l| are labelled by their energetic rank. The 
energy scale is in units of e. 



One implication of the method we use to generate our 
samples of stationary points is that higher-energy min- 
ima that were not used as starting points for transition 
state searches are unlikely to lie at the bottom of an 
MSB. These minima lie at the higher end of a pathway 
that was found in a search from a lower-energy minimum. 
Therefore, most of these minima are directly connected 
to a lower-energy minimum, and so cannot be at the 
bottom of an MSB. Occasionally, the pathways do not 
connect back to the starting minimum but to another 
unsearched minimum in the sample. Only in these rare 



instances can an unsearched minimum be at the bottom 
of an MSB in our sample. Therefore, although the dis- 
connectivity graphs based upon MSB's provide a more 
global picture of the PES than the graphs in the previ- 
ous section, they are still limited by the incompleteness 
of our samples. The graphs only provide a reliable pic- 
ture of the PES around the n sea rch lowest-energy minima, 
where n soar ch is the number of minima from which tran- 
sition state searches have been performed. 

For LJ13 there is only one MSB, reflecting its ideal 
funnel character and the remarkable degree of connec- 
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tivity. For LJig there are only a few MSB's, and these 
are all directly connected to the primary MSB, again re- 
flecting the single funnel character of the PES that we 
noted earlier. The double-funnel character of the LJ38 
PES is still apparent in the MSB disconnectivity graph, 
but now we get a better impression of the overall shape 
of the PES (Figure [j|a) . There is a wide, gently sloping 
funnel down from the higher-energy minima towards the 
low-energy icosahedral funnel, whereas the funnel down 
to the global minimum is much narrower. 2292 of the 
minima lie on monotonic sequences to the lowest-energy 
icosahedral minimum, whereas only 518 lie on sequences 



leading to the global minimum. Only 12 of the minima lie 
on sequences to both, showing that there is little overlap 
between the two MSB's. 

For LJ55 the MSB analysis leads to a remarkable sim- 
plification of the disconnectivity graph (Figure |b). The 
single-defect minima produce just one MSB, and the fine 
structure of the two-defect minima collapses onto the 
band of MSB's that branch off at -273e and -272e. The 
remaining branches are mostly three-defect minima with 
the three surface atoms close together. The graph clearly 
shows the single- funnel character of the PES. 
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FIG. 7. Disconnectivity graphs for the transformed energy surfaces of (a) LJ19 and (b) and (c) LJ38. In (a) and (b) 
only the branches corresponding to the lowest n ae arch steps are shown. In (c) only the branches corresponding to monotonic 
sequence basin bottoms are shown. The branches associated with the minima depicted in Figure |l| are labelled by their 
energetic rank. The energy scale is in units of e. 



C. Transforming the PES 

One approach to global optimization is to transform 
the PES to a form for which it is easier to find the 
global minimum.EJ One efficient-approach of this tyae is 
the Monte Carlo minimizationEil or basin-hoppingU al- 
gorithm. This unbiased method succeeded in finding the 
global|aninimum for all the clusters considered in this 
paperO In this approach the transformed potential en- 
ergy E c is given by 

£ C (X) =rnin{£ c (X)}. (2) 

Hence, the potential energy at any point in configura- 
tion space is assigned to that of the local minimum ob- 
tained by a minimization starting from that point, and 
the PES is mapped onto a set of interpenetrating stair- 
cases with steps corresponding to the basins of attraction 
surrounding a particular local minimum. Figure |^ illus- 
trates the transformation for a simple one-dimensional 



example. Note that the transformation does not change 
the identity of the global minimum, nor the relative ener- 
gies of any of the minima, but it does remove the downhill 
barriers between directly connected minima. This lat- 
ter change has a. significant effect on the dynamics and 

thermodynamics P^P 5 ! 

Due to our PES search strategy virtually all the un- 
searched higher-energy minima are directly connected to 
lower-energy minima. As the transformation removes the 
downhill barriers, in the disconnectivity graphs these un- 
searched minima would be directly connected to the stem 
associated with the superbasin which contains the rele- 
vant lower-energy minimum. This connectivity makes 
the higher-energy regions of the PES look funnel-like ir- 
respective of their actual character. Therefore, in the 
graphs we only depict branches which lead to the lowest 

^search Steps. 

The two examples shown in Figure |^ illustrate the ef- 
fects of the staircase transformation on the disconnectiv- 
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ity graphs. It removes the barriers to progress down a 
funnel, and the disconnectivity graph for LJ19 is therefore 
transformed into that for an ideal funnel. The long, dan- 
gling branches that are indicative of large barriers have 
disappeared, and so relaxation is now easier. However, 
barriers to interfunnel passage remain because the major 
component of such barriers usually arises from the high 
energy minima that the system has to pass through to 
go between funnels. (Figure ||a). For example, the trans- 
formation reduces the potential energy barriers for inter- 
funnel passage by only 0.68e for LJ 38 and 0.86e for LJ75. 
Therefore, the splitting of the LJ38 disconnectivity graph 
into two funnels is still clear, and perhaps even more ob- 
vious because many of the other barriers have been re- 
moved (Figure 0b). As for LJ19, the icosahedral region 
of configuration space now looks much more funnel-like. 




FIG. 8. A schematic diagram illustrating the effects of 
our potential energy transformation for a one-dimensional 
example. The solid line is the potential energy of the orig- 
inal surface and the dashed line is the transformed energy, 

The retention of the energy barriers for interfunnel pas- 
sage in LJ38 and LJ75 means that the global minima of 
these clusters-are still hard to find even on the trans- 
formed PES.c2l In fact, the success of the basin- hopping 
algorithm for these clusters lies in the changes-pto the 
thermodynamics caused by the transformation I 33 " 55 ! The 
thermodynamic transitions are broadened, thus provid- 
ing a temperature window where both the states at the 
top of the paths for interfunnel passage and the states at 
the bottom of the two funnels have a significant proba- 
bility of being occupied, making passage between funnels 
easier. 

On the transformed PES any pathway that is mono- 
tonically decreasing in energy (albeit a stepped rather 
than a smooth one) must end at the step corresponding 
to the bottom of an MSB. Therefore, an MSB on the 
transformed PES is analogous to a basin of attraction 
surrounding a minimum on the original PES. Figure ^ 
illustrates this point: on the transformed PES the two 
MSB's are like 'stepped' minima. However, the anal- 
ogy breaks down in one respect. A basin of attraction 
surrounding a minimum is defined as the set of points 
from which steepest-descent paths lead down to that min- 
imum. The use of the stcepcst-dcscent path ensures that 



each point in configuration space is mapped uniquely to 
a minimum even when a point is higher in energy than 
the lowest transition state connected to that minimum. 
However, on a step on the transformed PES the gradient 
is zero and so no steepest-descent direction is defined. 
Therefore, the mapping of a point in configuration space 
to a MSB bottom may not be unique if it lies above the 
divide between two MSB's, and so MSB's, unlike basins 
of attraction, can overlap. 

In Figure f?]c we show the disconnectivity graph for the 
transformed PES of LJ38 with only the branches leading 
to the bottoms of MSB's displayed. It has a similar struc- 
ture to Figure 0b but the density of branches is greatly 
reduced. 



IV. CONCLUSIONS 

In this paper we have shown that disconnectivity 
graphs can provide a valuable tool for visual represen- 
tation of a PES, and that they can provide clear physical 
insights into structure, dynamics and thermodynamics. 

The exponential increase in the number of minima 
with system size means that disconnectivity graphs with 
branches representing individual minima must, as the 
size increases, inevitably focus on low-energy regions 
of the PES which represent a decreasing proportion of 
the whole configuration space. By taking monotonic se- 
quences basins as the basic topographical unit, we have 
extended the cluster size for which a disconnectivity 
graph can represent the full sample of minima. However, 
these graphs still cannot provide a truly global picture of 
the PES because of the incompleteness of our samples of 
stationary points. 

We have also used disconnectivity graphs to probe the 
transformed PES that is searched by the Monte Carlo 
minimization or 'basin-hopping' global optimization al- 
gorithm. The transformation removes downhill barriers, 
making relaxation down a funnel much easier. For exam- 
ple, the algorithm can find the global minimum of LJ55, 
which has a single- funnel PES, in on average fewer &aa. 
150 steps when started from a random geometryoH 
However, the barriers between funnels remain, and so the 
success of basin- hopping for double- funnel PES's must be 
explained in terms of the different thermodynamics for 
the transformed landscape. l?JvZ3 
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